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Abstract. For real functions <J> and \P that are integrable and compactly 
supported, we prove the norm resolvent convergence, as e — > 0, of a family S e 
of one-dimensional Schrodinger operators on the line of the form 
d 2 

S E = -— - +ae- 2 *(e- 1 x) + /3e- 1 1'(e- 1 x). 
dx z 

The limit results are shape-dependent: regardless of the convergence of po- 
tentials in the sense of distributions the limit operator So exists and strongly 
depends on the pair (<t>, We show that it is impossible to assign just one 
self-adjoint operator to the pseudo-Hamiltonian — + ctS'(x) + 05(x), which 
is a symbolic notation only for a wide variety of quantum systems with quite 
different properties. 



1. Introduction 

The Schrodinger operators with singular potentials supported on a discrete set 
(such potentials are usually termed "point interactions") have attracted consid- 
erable attention both in the physical and mathematical literature from the early 
thirties of the last century. To understand the nature of quantum systems it ap- 
peared conceivable to analyze their general features about interactions with a range 
much smaller then the atomic size. Historically point interactions were introduced 
in quantum mechanics as limits of families of squeezed potentials. The quantum 
mechanical models that are based on the concept of zero range quantum inter- 
actions reveal an undoubted effectiveness whenever solvability together with non 
triviality is required. General references for this fascinating area are [1,5], which 
provide extensive documentation of pertinent material. 

In spite of all advantages of the solvable models, which are widely used in various 
applications to quantum physics, they give rise to many mathematical difficulties. 
One of the main difficulty of the analysis of zero-range interactions, compared to 
Schrodinger operators with short-range potentials, is that the Schrodinger opera- 
tors with singular potentials are often only formal differential expressions, and for 
the corresponding differential equations no solution exists even in the sense of dis- 
tributions. In 1961 Berezin and Faddeev [8] suggested how such formal Schrodinger 
operators can be constructed as mathematically well-defined objects, and for the 
first time a formal Hamiltonian was written as a self-adjoint operator derived by the 
theory of self-adjoint extensions of symmetric operators (see [20] for more details). 
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There exists a large body of results on this subject. It is impossible to refer to 
all relevant papers, and I confine myself to a brief overview of the one-dimensional 
case. In the recent years, many results were obtained for the point interactions 
based on the theories of self-adjoint extensions of symmetric operators, singular 
quadratic forms, boundary triples and almost solvable extensions [4,25-28], and 
once more we refer to [1, 5] for a exhaustive list of references therein. Another 
way to define a formal Schrodinger operator with a distributional potential v £ 
2?'(R) is to approximate it by Schrodinger operators with more smooth potentials 
v £ obtained by a suitable regularization as well as to use the concept of quasi- 
derivatives [6,13,15,21,22,32-35]. In [14,16,29,31] the solvable models for ID 
Schrodinger operators were based on specific products of 5^ and discontinuous 
functions, where 6 is the Dirac delta function. It is worth to note that the difficulty 
of dealing with the multiplication in V may also be overcome by using the new 
algebras of generalized functions [7] . 

It is common knowledge that all nontrivial point interactions at a point x can 
be described by the coupling conditions 

flp(x + 0)\ _ _ fl/)(x - 0)\ r _ c i V ( c ll c 12\ /"i i \ 

y(x+o))- c y(x-o))> c ~ e [ C21 c 22 )> 

where ip 6 [— §-,§], cti G K and c\\c 2 2 — c\ic%\ — 1< For the physically based 
classification of these interactions we refer the reader to the recent preprint [9], 
where in particular one singles out the four most important cases: (^-potentials, 
S -interactions, ^'-potentials and (^-magnetic potentials. For a quantum system 
described by the Schrodinger operator with a smooth enough potential localized in 
a neighbourhood of x one can often assign a point interaction with some matrix in 
(1.1) so that the corresponding zero-range model governs the quantum dynamics 
of the true interaction with admissible fidelity, especially for low-energy particles. 
However the connection between real short-range and idealized point interactions is 
very complex and ambiguously determined. This is certainly the reason why there 
is a number of papers on this subject occasionally even with conflicting conclusions. 

As for the (5-potential, any smooth approximation of j38[x) leads to the same 
solvable model given by conditions (1.1) with the matrix 

c-(j;). M 

Thus the result is shape-independent, i.e., it is not sensitive to regularization. In this 
case the limiting argument admits a straightforward interpretation. The nonzero 
off-diagonal element of C implicitly involves the integral of the approximating po- 
tential, and hence a slow particle on the line "feels" only the average value of a 
localized potential. 

The situation changes if we turn to the (^'-potential. The usual regularization 
of S'(x) is a sequence e~ 2 v(e~ 1 x) with a zero-mean function v € Cq°(R). It was 
shown in [18] that for almost all functions v the best zero-range approximation to 
the Hamiltonian H e = — j^i + e~ 2 v(s~ 1 x) is the free Hamiltonian —^2 subject 
to the split boundary conditions 0) = i/j(+0) = 0. These conditions define a 
non-transparent interaction at the origin. However, there exist so-called resonant 
potentials v (see below for the precise definition) for which the limit behaviour of 
quantum system can be characterized by the nontrivial point interaction with the 
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coupling matrix 

where = 9{v) is a spectral characteristics of the potential v. Incidentally, it is 
of interest that for any shape v there exists a countable set of so-called resonant 
coupling constants for which the spreading potentials a^v are resonant. These 
results were recently extended to potentials v of the Faddeev-Marchenko class [19] 
and generalized to the case of quantum graphs [30]. 

Therefore the results on ^'-potentials become shape-dependent: depending on v 
the Hamiltonians H e regularize different kinds of point interactions, nevertheless 
all of them involve the 5'-like potentials. Hence, it is impossible to assign just one 
self-adjoint operator to the pseudo-Hamiltonian — ^ +aS'(x), which is a symbolic 
notation for a wide variety of quantum systems with quite different qualitative and 
quantitative characteristics. 

It has been known for a very long time that the (^'-potential defined through 
the regularization e~ 2 v(e~ 1 x) is opaque acting as a perfect wall, see widely cited 
Seba's paper [35] of 1986. However, such a conclusion is in contradiction with 
the analysis of H e with piece-wise constant potentials v performed recently by 
Zolotaryuk a.o. [11,36-38], where the resonances in the transmission probability 
for the scattering problem are established. In [17] a similar resonance phenomenon 
is also obtained in the asymptotics of eigenvalues for the Schrodinger operators 
perturbed by <5'-like potentials. The authors of [2, 10] faced the question on the 
convergence of H e in approximation of a smooth planar quantum waveguide with 
a quantum graph. Under the assumption that the mean value of v is different from 
zero, they also singled out the set of resonant potentials v producing a non-trivial 
limit of H £ in the norm resolvent sense. The situation with these controversial 
results was clarified in [18]. Curiously enough, P. Seba was the first to discover in 
1985 [34] the resonant potentials for a similar family of the Dirichlet Schrodinger 
operators on the half-line producing in the limit the Robin boundary condition. 

This paper can be viewed as a natural continuation of the recent work [17-19] on 
the Schrodinger operators with (5'-likc potentials to the case in which the potentials 
are a smooth enough regularization of the distribution aS'(x) + f}6(x). Clearly it is 
to be expected that the limit results concerning such families of squeezed potentials 
will be also shape-dependent. 

Notation. Throughout the paper, W^oj) stands for the Sobolev space of functions 
defined on a set oj C M that belong to Lzfa) together with their derivatives up to 

1/2 

order I. The norm in Wf (w) is given by ||/||wf ( fi) = (||/"||£ 2(n ) + ll/ll| 2 (n)) . 
where ||/||t 2 m) is the usual L 2 -norm. We shall write ||/|| instead of ||/||l 2 (r). 

2. Statement of Problem and Main Result 

Let us consider the formal Hamiltonian 

d 2 

H=-—+a6'{x)+l36(x), x € R, 
dx z 

where 5 is the Dirac delta function. In accordance with the classic theory of dis- 
tributions we have (H(f>){x) = -<f>"(x) + a<j>(Q)6'(x) + 080(0) - a<f>'(0))6(x) for any 
4> € C 1 (M). However, there exist no solutions in V (R) to the equation Hcf> = \<f> 
for a nonzero a, except for a trivial one. The reason for this at first sight surprising 



4 



YURIY GOLOVATY 



fact lies in the theory of distributions: X>'(R.) is not an algebra with respect to the 
"pointwise" multiplication. 

Instead of H, we consider the family of Schrodinger operators 

s -=-^ + ;H?) + i*(?)- < 2 -d 

with integrable potentials $ and iff of compact supports. Here £ is a small positive 
parameter, and the coupling constants a and /? are assumed to be real. One of the 
questions of our primary interest in this paper is the behaviour of S s as e tends to 
zero. The motivation for this question stems from the fact that the potentials 

V e = as- 2 ®^- 1 ■ ) + /fe- 1 *^ 1 • ) 

approximate the pseudopotential aS' + j36 under some assumptions on $ and iff. 
Indeed, if the following conditions hold 

<f>ds = 0, [s$ds = -l and f^ds = l, (2.2) 
Jr Jr Jr 

then ae~ 2 &(e~ 1 x) + fie~ l if!{e~ l x) — > aS'(x) + (35(x) as e — > in the sense of 
distributions. In this case, we call <E> the shape of a <5'-like sequence and W the 
shape of a 5-like one. 

Notwithstanding the title of paper, all results presented here concern the poten- 
tials V E with arbitrary $ and iff of compact support, and the aS' + /35-like potentials 
are only a partial case in our considerations. Note that if the first condition in (2.2) 
is not fulfilled, then the potentials V s do not converge even in the distributional 
sense. However, surprisingly enough, regardless of the convergence of V E the limit 
of S £ exists in the norm resolvent sense (i.e., in the sense of uniform convergence 
of resolvents). 

Definition 2.1 ([24]). We say that the Schrodinger operator —j^i + q in L2(K) 
possess the half-bound state (or zero-energy resonance^ provided there exists a 
solution u to the equation —u" +qu = inM. that is bounded on the whole line, i.e. 
u G L°°(R). The potential q is then called resonant. 

Such a solution u is unique up to a scalar factor and has nonzero limits u(±oo). 
Our main result reads as follows. 

Main result. Let $ and iff be integrable and bounded real functions of compact 
support. Then the operator family S e given by (2.1) converges as e — > in the norm 
resolvent sense. 

If the potential a$ is resonant with a half-bound state u a> and u~ = u a (±oo), 
then the limit operator So is a perturbation of the free Schrodinger operator defined 
by S (j) = -4>" on functions (f> in W 2 2 (R \ {0}) obeying the boundary conditions at 
the origin 

</>(+0) - a <t>(-0) = 0, (f>'(+0) - 9~ V (-0) = [3 x a <f>(-0). (2.3) 
The parameters a and >t a are specified by the potentials $ and ^ : 

u + If 

6 a = -§h Xa = / f>U 2 a dt. 

U a U a U a JR 

Otherwise, in the non-resonant case, the limit S$ is equal to the direct sum 
5_ ® S+ of the Dirichlet half-line Schrodinger operators S±. 
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The result is proved in Theorems 4.1 and 5.1 below. In addition, convergence of 
the scattering data for S e to the ones for So is established in Theorem 6.1. 
In the resonant case, the point interaction generated by the coupling matrix 





C7($,*) = 



13 K a 



3-1 




in (1.1) may be regarded as a first approximation to the real interaction governed 
by the Hamiltonian S £ . The explicit relation between 9 a , >c a and the potentials 
^> makes it possible to carry out a quantitative analysis of the quantum system, for 
instance, to compute approximate values of the scattering data for given $ and <!/. 

It is appropriate to mention here that in [3, 29] the pseudo-potential ad' + [35 
was interpreted as a point interaction with the matrix 

C= 

\(2-«) 2 

and some split boundary conditions were associated with the singular values a = 
±2. The solvable model was derived from the assumption that the following product 
formulae 

w(aO*(a;) - Mo S(x), v{x)S'(x) = {v} S'(x) - {v'} S(x) 

hold, where {/}o = |(/(— 0) + /(+0)) is the mean value of a discontinuous function 
/ at x = 0. The spectrum and scattering properties of this model were described 
in [14,16]. 

3. Resonant Sets and Maps 

Since the potential V e has compact support shrinking to the origin, there is no 
loss of generality in supposing that the supports both of <1> and "3/ are subsets of the 
interval / = [—1, 1]. Denote by V the class of real integrable and bounded functions 
of compact support contained in /. 

Definition 3.1. The resonant set A$ of potential $ £ V is the set of all real value 
a for which the operator + ct^ in L 2 (IR) possesses a half-bound state. 

Suppose that a potential q g V is resonant, i.e. q possesses a half-bound state u. 
Then u, as a solution to the equation —u" + qu = 0, is constant for \s\ > 1, because 
g is a zero function outside /. Moreover, the restriction of u to / is a nontrivial 
solution of the problem — u" + qu = 0, s € (—1,1) and u'{— 1) = 0, u'(l) = 0. 
Hence, the potential q is resonant if and only if zero is an eigenvalue of the operator 
N = —-^2 + q in L2(I) subject to the Neumann boundary conditions at s = ±1. 

Consequently the resonant set A$ coincides with the set of eigenvalues of the 
problem 

-u" + a$u = 0, se(-l,l), it'(-l) = 0, u'(l) = (3.1) 

with respect to the spectral parameter a. In the case of a positive $ it is clear 
that A$ is a countable subset of K+ without finite accumulation points, and all 
eigenvalues of (3.1) are simple. Otherwise, (3.1) is a problem with indefinite weight 
function [12]. 

Assume that $ has only isolated turning points in /. This case was considered in 
[17]. Then problem (3.1) can be associated with an operator in an appropriate Krein 
space. Let /C$ be the weight L 2 -space with the scalar product (/, g) — J_ 1 fg |$|ds 
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and the indefinite inner product [/, g] = (J/, g), where J f = sgn $ • /. The operator 
J is called the fundamental symmetry. We can introduce in /C<j> the operator 

T = ~$^)S' ^T={gelC^.geWi(I), 

$-y'e£*, g'(-i) = 0, g'(i) = 0} 

that is J-self-adjoint and J-nonnegative. The spectrum of T is real and discrete, 
and has the two accumulation points — oo and +00. All nonzero eigenvalues are 
simple, and a — is semi-simple, generically. The reader is referred to [23] for 
the details of the theory. Obviously, A$ = — <r(T). Hence, the resonant set A$ is 
discrete and unbounded in both directions. 

Now suppose that the support of $ is a disconnected subset of I. For the sake 
of simplicity, assume that it has only one gap: supp$ = [— l,si] U [s2,l] and 
S\ < s 2 . Each solution u of (3.1) is then a linear function on [s\,S2]- Therefore 
u'(s 2 ) = u'(si) and u(s 2 ) — u(s±) = lu'(si), where I = S2 — s± is the length of gap. 
Let us move the interval [— l,s\] up to [s2,l], identify the points s± and S2, and 
thereafter rewrite problem (3.1) as 

W + qTh = 0, x G [I - l,s 2 ) U (s 2 ,l), v'(l-l)=0, v'(l)=0, 

v'(s 2 +0)=v'(s 2 -0), v{s 2 +0)-v(s 2 -0)=lv'{s 2 ). '' >J 

The new "glued" potential T coincides with <f> on [s 2 , 1], T(s) = $(s — I) for 
s £ [I - l,s 2 ); hence T has a unique turning point s = s 2 . The relation between 
solutions of problems (3.1) and (3.2) is obviously given by 

{v(s + I) for s G [-Mi), 

v'(s 2 )(s - si) + u(s 2 - 0) for s G [si,s 2 ], 
u(s) for s G (s 2 , 1). 

As in the previous case, we can now construct an J-self-adjoint and J-nonnegative 
operator in tCy associated with problem (3.2) and derive the same properties of the 
resonant set A$. The similar considerations can be applied to $ with several gaps 
on its support. 

In conclusion of this section, we introduce two characteristics of the potentials 
$ and which will turn out to be important for us later. Let us denote by u a the 
half-bound state that corresponds to resonant potential a<I>. Clearly, u a (±oo) — 
it a (±l). Let 6 be the map of A$ to K such that 

9(a) = — — - (3.3) 
u a (-l) 

for all a £ A$. Next, let the map x is given by 

x{a) = 1 / ^u 2 a da, a G A«. (3.4) 

The value x(a) describes the interaction of potentials $ and * at the resonant 
a. Since the half-bound state is unique up to a scalar factor, both maps are well 
defined. 



Definition 3.2. We call 6: A$ — > K the resonant map 0/$, and x: A$ — > K i/ie 
intercoupling map for a pair of potentials $ and . 
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4. Convergence of S s in Resonance Case 

In this section, we analyze more difficult resonant case where a is a point of the 
resonant set A$. Here and subsequently, a = 9(a) and x a = x(a), where 9 and 
x are the resonant and intercoupling maps for a given pair of $ and ^> . Let u a be 
the eigenfunction of (3.1) corresponding to a £ A$ such that u a (— 1) = 1. From 
(3.3) and (3.4) it follows that 9 a = u a (l) and 

x a = e- 1 J *u 2 a dt. 

Denote by S(fj,, v) the free Schrodinger operator on the line acting via S(fi, v)<p = 
—<j)" on the domain 

dam S(ji,v) = {<p€ W|(K\{0}): (f>(+0) = fuf>(-0), <j>'(+0) = /i'V (-O)+z/0(-O)}. 
For each real v and fi ^ the operator S(p, v) is self-adjoint. 

Theorem 4.1. Assume that G V and a belongs to the resonant set A$. Then 
the operator family S s defined by (2.1) converges to S(9 a , (3x a ) as s — > in the 
norm resolvent sense. 

We have divided the proof into a sequence of propositions. Fix an arbitrary 
/ G L20R-) and ( e C with Im£ ^ 0. The basic idea of the proof is to construct a 
fair approximation to the function y £ — (S s — Q f, uniformly for / in bounded 
subsets of L20R). 

In the sequel, letters Cj, Cj and bj denote various positive constants independent 
of e and /, whose values might be different in different proofs, and ||/|| stands for 
the i2(K)-norm of a function /. For abbreviation, we let So stand for S(9 a , f3x a ). 

Set y = (So — C) 1 f ■ We show that y is a very satisfactory approximation to y e 
for \x\ > e. The problem of choosing a close approximation to y £ on the support of 
potential V e is more subtle. Denote by v e the solution of the Cauchy problem 

(-i% + a$(s)v e =ef(es)-Py(-e)V(s)u a (s), s € (-1,1), 

\ v e (-l) = 0, v >(-l)=y>(- £ ). [ ■ > 

Clearly, we have v £ S W|(— 1, 1) for any e > 0. The next proposition establishes 
some asymptotic properties of v e as e — > 0. 

Proposition 4.2. The following holds for any f G L 2 (M.) and e G (0, 1) 

IMIw?(-u) < CiWfl K0-) - y'(+o)l < <V /3 II/||. 

Proof. We first observe that (So — C) -1 is a bounded operator from L 2 (M.) to the 
domain of So equipped with the graph norm. The latter space is a subspace of 
Wf (M \ {0}), and therefore 

IMIw?(h\{o}) < cxll/ll, lbllci(K\{0}) < c 2 ||/||, (4.2) 
since Wf (M \ {0}) C C 1 (R\ {0}) by the Sobolev embedding theorem. Hence, 

IKIkf(-i,i) < c 3 (\y(-e)\ + \y'(-e)\ +e||/(e-)IU a (-i,i)) < C ^\fl 
because there exists a constant C4 such that 

||/( £ -)IU 2 (-i,i)<C4 £ - 1/2 ||/||. (4.3) 
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Next, multiplying equation (4.1) by the eigenfunction u a and integrating by parts 
yield 

<(1) = 9- l y'(-e) + f3x a y(-e) - e^ 1 J ^ f(es) u a (s) ds. (4.4) 

Recall that the function y satisfies the condition y'(+0) = 0~ 1 y'(— 0) + {3x a y(—0). 
Subtracting this equality from (4.4) we finally obtain 

|<(1) - y'(+0)\ < IB" 1 ] \y'(-e) - y'(-0)\ + |/3| |x Q | \y(-s) y(-0)\ 

+ ||/(e-)lk(-i,i)IKIU 3 (-i,i) < C2S 1/2 \\f\\ 

in view of (4.3) and the following estimates 



LW(± £ )- y W(±0)|< 



±e 



\y {k+1 \x)\dx 



(4.5) 

holding for k = 0, 1. □ 

Let us introduce the function w e such that w £ {x) = y{x) for |x| > e and w e (x) = 
y(—e)u a (x/e) + ev E {x/e) for \x\ < e. By construction, w e belongs to W$ (K \ {s}). 
Although w £ is in general discontinuous at the point x = e, its jump and the jump of 
its first derivative at this point are small. Indeed, [u> e ]a; =£ = y{s) — 6 a y{— e) — sv e (l) 
and [w' e ] x=e = y'(e) — v' s (l), where [/i]a; =a = h(a+0)— h(a— 0) is a jump of a function 
h at x = a. Therefore, taking into account Proposition 4.2, estimates (4.5) and the 
equality y{+0) = 9 a y(—0), we see that the jumps can be bounded as 

|[«4*=«| < \y(e) - V(+0)\ + \6 a \\y(-e) - y(-0)\ + s\\v E \\ wi{ _ 1A) < 6i£ 1/2 ||/||, 
|KW| < \y'(e) - y'(+0)\ + \y'(+0) - <(1)| < he^Wfl 

(4.6) 

Let us introduce functions r/o and r]\ that are smooth outside the origin, have 
compact supports contained in [0,1], and have the prescribed jumps [rjo} x =o = 1, 
Wo\x=o = and fri]^ = 0, [?7j] x= o = 1 (sec Fig. 1). Set z e {x) = [w e ] x=E r] {x - 
e) + [w' e ] x=E rji (x — e); then in view of (4.6) 

maxIzW^l^e 1 / 2 !!/!! (4.7) 

x>s 

for some 63 and k = 0, 1, 2. Additionally, z £ = on (—00, e). 
Clearly, the function 

- M = f 2/0*0 -z e {x) ii\x\>e, 

Ve[X> \y(- e )u a {2)+ev a (1) i£\x\<e { ' ' 

is continuous on R along with its derivative and belongs to dom S £ . 




Figure 1. Plots of the functions with the prescribed jumps at the origin 
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Proposition 4.3. Fix ( e C with Im£ ^ 0. Then the estimates 

life -fell <Cie 1/2 |||/H, \%-v\\<C^ 2 \\f\\ (4.9) 
hold for each f G and e > 0, where y e — (S £ — C) -1 / an d V = (So — C) / • 

Proof. It is convenient now to rewrite the approximation y £ in the form 

y e (x) = (l-x s (x))y(x)+y(-e)u a (^) + sv e (^)-z s (x), 

where \e is the characteristic function of [— e,e], and u a and v e are extended by 
zero to the whole line. Recalling the definition of y, u a and v e , we deduce: 

(S e - Qy e (x) = (-£ - 0(y(x) - z e (x)) = f(x) + z'J(x) + (z e (x) 
for |x| > e, and 

(S s -Qy e (x) = (-£, + ae-^e-'x) + ^-^(e^x) ~ c) (l/(-e)«a(f ) + £»«(§)) 

= e- 2 y(- £ ){-< + a*(f )u a } + £- x {-< + «*(f H + /%(-e)*(f 
+ /?*(§ )„ e (§ ) - Cy e (x) = /(a,) + (/?*(§) - C)fe(f) - Cw(-e)«a(f ) 
for |x| < e. Therefore (S e — C)fe = / + r e , where 



r 6 (ac) 



f4'(a:) + C^(a;) if |x| > e, 

\ (/3*(f ) - <K(f ) ~ Cy(-eK(f ) if M < e. 



Hence y e — y s = (S £ — C) r e) and from this we conclude 



|y e -£ e ||<||(S £ -Cril|re||<|ImCr ] 



r F 



We can now employ Proposition 4.2 and estimates (4.2), (4.7) to derive the bound 

\\r e \\ < ci||z"+Cz e || L2(e4+e) +c 2 ||w e (e _1 • )\\ L2 (- e , e )+c 3 \y(-e)\ ||tt a (e _1 • )||i 3 (_ e , e ) 
< C4max(|z e | + |4'|)+c 5 e 1/2 (||^|| L2( _ l!l) + ||y|| c(R \ {0}) ||u Q || L2( _ lil) ) < c 6 e 1/2 ||/||. 

This proves the first inequality in (4.9). Similarly, 

\\ye-y\\ = \\y(-£)u a (e~ l -)+ev e (e- l -)-z s -Xey\\ < c 7 e 1/2 |?;(-e)|||Ma|U 2 (-i,i) 
+ c 8 e 3/2 ||w e || L2( _ 14) +c 9 max| 2£ | +c 10 ||j/|| c(R \ {0}) ||x £ || < cii£ 1/2 ||/||, 

and so finish the proof. □ 

Proof of Theorem 4-1- For each / e L 2 (M.) and ( £ C with Im£ ^ we can 
construct the approximation y e to y e = (S e — C) 1 f given by (4.8). As above, 
set y = (So — C) -1 /- Applying Proposition 4.3, we discover 

\\(S e - c)- 1 / - (So crVll = life - y\\ < life - fell + life - y\\ < Ce l / 2 \\fl 

which establishes the norm resolvent convergence of S e to the operator S(9 a , f3x a ). 

□ 
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5. Convergence of S e in Non-Resonance Case 

Now we study the non-resonant case when a does not belong to the resonant set 
A$ . Denote by Sq the direct sum S— © S+ of the unperturbed half-line Schrbdinger 
operators S± = —d 2 /dx 2 on R± subject to the Dirichlet boundary condition at 
x = 0. Hence domSo = {y € W${R \ {0}) : y(-0) = y(+0) = 0}. 

Theorem 5.1. If a ^ A$, i/ien t/ie operator family S e given by (2.1) converges to 
the direct sum S- © S+, as e — > 0, in the norm resolvent sense. 

Proof. Exactly the same considerations, as in the previous section, apply here, with 
one important difference: the function y = (So — £)/ is small in a neighbourhood 
of the origin, since y(0) = 0; and we have to change the approximation to y £ — 
(S e — C) -1 / on the support of V e . We set y e (x) — y(x) — z e (x) for \x\ > e and 
y e (x) = ev e (~) for | a; | < e, where v e is a solution to the boundary value problem 

- < + a$(s)v s = ef(es), s e (-1, 1), t/ e (-l) = ^(-e), t/ e (l) = y'(e), (5.1) 

and z e (a;) = (e« s (— 1) — y(—e)J ?yo(— x — e) + (j/(e) — eu e (l)) r] (x — e). Note that 
the solution v £ of (5.1) exists, because in the non-resonance case the number a is 
not an eigenvalue of the corresponding homogeneous problem. Moreover, 

Ikll^K-i.i) < c x (\y'{-e)\ + \y'(e)\ + - )IU a (-i,i)) < c 2 ||/||, 
due to (4.3) and the apparent estimates ||y||ci(B\{o}) < c 3||2/|Ivk 2 2 (e\{o}) < c 4||/ll- 

Recalling that y(0) = 0, we deduce from (4.5) that |y(±e)| < c 5 e 1/2 \\f\\. There- 
fore the corrector z e can be bounded as 

maxl^^l^cee 1 / 2 !!/!! 

\x\>e 

for k = 0,1,2, since \y{±e) - ev e (±l)\ < |y(±e)|+e|v s (±l)| < ere 1 / 2 1| /|| . 
Next, an easy computation shows that (S e — Qy e = / + where 

f<(or) + Cz £ (x) ifH>e, 
reW \(Wf)-<K(f) ifW<e. 

As in Proposition 4.3, from the estimates above we can derive ||r e || < cg£ 1 //2 1 1 ^ 1 1 - 
Hence, 

|| (5 S - O^f - (So - O^fW <\\Ve- tf e || + \\Ve- V\\ 

< \\(S e ~ CrV|| + Wev^e- 1 -)-z £ - xey\\ < Ce l / 2 \\f\\. : 
which establishes the norm resolvent convergence of S e to the operator £_ ®S+. □ 

6. Scattering on a5' + /3<5-Like Potentials 

6.1. Scattering problem for Sq. First, let us discuss stationary scattering asso- 
ciated with the Hamiltonians S(9 a , 0x a ) and — ^2, which corresponds to the reso- 
nant case. Consider the incoming monochromatic wave e lkx with k > coming from 
the left. Then the corresponding wave function has the form ip(x, k) = e lkx +R e~ ikx 
for x < and ip(x, k) = T e lkx for x > 0. Here R and T are respectively the re- 
flection and transmission coefficients. The matching conditions (2.3) at the origin 
clearly yield 

(T\_(9 a \ ( 1 + R \ 

\ikT) \pK a e- 1 ) [ik(i-R)J ■ 
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Then one obtains the reflection and transmission coefficients from the left, cf. [16, 
Eq 23]: 



a) 



fix a 



T(k,a) 



2ik 



fix a 



a G Ad 



(6-1) 

In the non-resonant case the scattering problem is trivial. The split condition 
y(0) = leads to the equalities R(k, a) = — 1 and T(k, a) = 0. 

6.2. Convergence of the scattering data. Next we investigate the stationary 

,2 

scattering for S e (a, (3, <&, \I/) and — , and prove that the scattering data converge 
as e — > to the scattering data for So obtained above. We look for the positive- 
energy solution to the equation — tp" + V E ip = k 2 ip given in the form 



e'"" + R E e~ lKx for x < -s, 

A £ u e (e~ 1 x, a) + B e v e (e~ 1 x, a) for |x| < e, 
T F e ihx for x > e. 



Here u s = u e (s, a) and v e = v e (s, a) are solutions of the equation 
- w" + a$(s)w + Ps^(s)w = e 2 k 2 w 1 s G (-1, 1) 
subject to the initial conditions 

u e {— 1, a) = 1, u' e (— 1, a) = and V £ (— l,a)=0, v' e (—l,a) 



1 



(6.2) 



(6.3) 



respectively. The coefficients R E , A e , B £ , and T e can be found from the linear 
system 





1 





° \ 


(Re\ 




( Z- {ek \ 


ieke 1£ 





1 





A e 




\eke 1£ 





u e (l,a) 


v e (l,a) 


-e lek 









V o 


<(l,a) 


<(l,a) 


-ieke lek J 






V o J 



obtained by matching the solution and its first derivative at the points x = ±e. By 
Cramer's rule, we can derive 



R s {k, a) 
T £ {k,a) ■■ 



,-2iefe 



-2iefe 



4(1, a) - iek(ue(l, a) - v' e (l, a)) + e 2 k 2 v e (l, a) 
u' £ (l, a) — \Ek(u e {X, a) + v' e (l, a)) — e 2 k 2 v E (l, a) ' 
2iefc 

u' e (l, a) — iek(u E (l, a) + v' e (l, a)) — e 2 k 2 v E (l, a) 



(6.4) 



Here we use the identity u E (l,a)v' £ (l,a) — u' £ (l,a)v E (l,a) — 1 that follows from 
(6.3) and the constancy in £ and e of the Wronskian of u E and v E . 

Theorem 6.1. For each k > and a G M the scattering data R E (k, a) and T E (k, a) 
converge respectively to i?(fc, a) and T(k, a) as e — > 0, where the limit values are 
given by (6.1) in the resonant case, and R(k,a) = —1, T(k,a) — otherwise. 

Proof. From the smooth dependence of a solution to the Cauchy problem on para- 
meters, we see that u E and v E converge in C 1 (— 1,1) to the solutions u and v 
respectively of the equation — w" + a$w = subject to the initial conditions 

u(— 1, a) = l, u'(— l,a) = and v(—l,a) = 0, v'(—l,a) = l. 
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i LA) 








a -' 





a, 


a 2 


a. 



Figure 2. A plot of transmission probability \T e (k, a)\ 2 as func- 
tion of a. 



The non-resonant case. Since a is not a eigenvalue of problem (3.1), we conclude 
that u'(l, a) is different from 0. From (6.4), it immediately follows that R e (k, a) = 
-1 + 0(e) and T £ (k, a) = 0(e) as e -> 0. 

TVie resonant case. If a is resonant, then u is an eigenfunction of (3.1). Therefore 
u'(l, a) = and u(l, a) = Throughout the proof, 8 and x denote the resonant 
and intercoupling maps for a pair ($, ^1/). Let us substitute functions u e and v e into 
(6.2) alternately. Multiplying the derived identities by u and integrating by parts 
yield 

9 a u' s (l,a)= e/3 ^u £ ud£ + e 2 k 2 u e udt;, 

,i /i (6-5) 

6 a v' s (l,a) = l + e@ J ^v E ud£_ + e 2 k 2 J v E ud£,. 

Therefore u' E (l, a) = e(3x a + o(e) and v' e (l, a) = O^ 1 + 0(e) as e — > 0. Combining 
then these asymptotic formulae and (6.4), we have 

R e {k, a )= ' lk{e \~ 6a)+PXa + (1), T e (fc,a) = ^ +o(l), 

(6.6) 

as e — > 0, which completes the proof. □ 



6.3. Resonances in the transmission probability. It follows from the theorem 
above that the probability of transmission across the barrier V £ = ae^ 2 ^(e^ 1 •) + 
fie~ Xi $ (e~ x ■) is negligibly small for a £ A$. However, for the resonant coupling 
constants a this probability remains non-zero as s — > 0, resulting in the existence 
of non-separated states. The resonances of transmission, as shown in Fig. 2, are 
sharp like one-point spikes, which spread for non-zero values as e — > 0. The limit 
resonant values can be represented via the maps 9 and x: 

Ak 2 

\T(k,a)\ 2 = = , a€A $ . (6.7) 

1 V 71 k 2 ^ 1 +6 a ) 2 + (3 2 x 2 a ' * V ' 

Several special cases are of interest. 

6.3.1. 6 -potentials. Set a = 0, then V e — j3e~ 1 '$(e~ 1 -). Assume also that J R ^> dt = 
1. Consequently V e converges to f3S(x) in the sense of distributions. Note that 
the value a — belongs to A$ for each $. The corresponding eigenfunction uq 
of (3.1) is constant on [—1,1], and both numbers 9q and are equal to 1 (see 
(3.3), (3.4)). Hence, the limit operator S(l,j3) describes the point interaction at 
the origin with the coupling matrix given by (1.2). In addition, from (6.7) we 
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have |T(fc,0)| 2 = 4fc 2+a2 , which is common knowledge of the scattering by the 
/3(5(a;)-potential. 

6.3.2. 5' -like potentials. Now we set /3 = 0. Then V £ = a£ _2 $(e _1 -). If J R dt = 0, 
then 

and so in particular the limit can be zero. Otherwise, the family of potentials V e 
diverges in the sense of distributions. Regardless of the weak convergence of V E the 
transmission probability |T e (£;,a)| 2 for each a € A$ converges to the value 

\T(k,a)\ 2 = =— (6.8) 

that does not depend on energy of particles (cf. [18, Sec. 5]). The limit operator 
S(a, 0) corresponds to the point interaction at the origin with the coupling matrix 
given by (1.3), 9 — 9 a . As opposed to the case of the ^-potential, this result is 
shape-dependent via the resonant map 9. 

Our analysis of the exactly solvable models with piecewise-constant potentials 
and the computer simulation of more complicated models may suggest that if 
e _2 $(e -1 -) — > 5'{x) as e — > 0, then \8 a \ — > +oo as a — > +oo and \8 a \ — > as 
a — > — oo. Therefore the transmission probability for the <5'-like potentials given by 
(6.8) is very small for large \a\. As we will see in the next special case, it is not 
true in general. 

It is noteworthy that the operator S(a, 0) can also appear as a solvable model for 
the Schrodinger operator with the potential V £ when /3 is different from 0. Indeed, 
for some resonant values of a it is possible for potentials $ and 4" to be uncoupled, 
i.e., K a = 0. If, for instance, a potential \& is sign-changing, then it is possible for 
the integral J _ x 'Sit 2 ds to be zero. 

6.3.3. Potentials with total transparency at resonances. An interesting case occurs 
when for a pair of potentials $ and the intercoupling map x is identically zero 
and the resonant map 9 is unimodular, \9\ = 1. Then the marginal transmission 
probability across the potential V e is given by 

1 if a G A$, 
otherwise. 

Hence, either the potential is asymptotically opaque for particles or else asymptot- 
ically totally transparent at resonances. For instance, this kind of case arises when 
$ is an even function, whereas \& is an odd one. Then each eigenfunction u a of 
(3.1) is either odd or even. In any case, \9 a \ = 1 and x a = 9~ x f_ t 'Ju 2 ds = 0, 
since the square of u a is an even function. 
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